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We investigate one dimensional tight binding model in the presence of a correlated binary disorder. 

The disorder is due to the interaction of particles with heavy immobile other species. Off-diagonal 
disorder is created by means of a fast periodic modulation of interspecies interaction. The method 
based on transfer matrix techniques allows us to calculate the energies of extended modes in the 
correlated binary disorder. We focus on A-mer correlations and regain known results for the case of 
purely diagonal disorder. For off-diagonal disorder we find resonant energies. We discuss ambiguous 
properties of those states and compare analytical results with numerical calculations. Separately we 
describe a special case of the dual random dimer model. 


I. INTRODUCTION 

In one dimensional disordered systems particles can¬ 
not freely propagate. The particle undergoes a series of 
scattering processes and, as a result of the destructive in¬ 
terference, its density profile decays exponentially. This 
phenomenon is called the Anderson localization (AL) 
[1, 2]. AL is a single particle effect and therefore can¬ 
not be observed directly in solids (because of the strong 
electron-electron and electron phonon interactions). 

Nevertheless, nowadays many different tight-binding 
models can be easily implemented and investigated with 
ultracold atoms in optical lattices. Ultracold atomic sys¬ 
tems in optical lattices guarantee a high degree of con¬ 
trollability [3-5], allowing to change experimentally pa¬ 
rameters of the model, such as hopping amplitudes using 
periodic lattice modulations [6-8], or interparticle inter¬ 
action strength using the magnetic Feshbach resonance 
[9[. In particular, the interactions can be switched off 
completely, creating perfect conditions for the investiga¬ 
tion of AL. Indeed, AL of matter waves has been already 
observed in the laboratory [10, 11] (as well as a similar 
Aubry-Andre localization [12, 13[). 

It is very convenient to study disorder phenomena with 
optical lattices. In particular, the diagonal disorder in ul¬ 
tracold systems can be very easily implemented by ma¬ 
nipulating on-site energies, e.g. with incommensurate 
lattices [13] or interactions with other (heavy) species 
of particles [14-18]. The later method allows to intro¬ 
duce various types of correlations in the disorder distri¬ 
bution, e.g. dual random dimer models (DRDM), where 
two heavy atom never come in pairs) or A^-mers, where 
heavy particles always occupy a series of neighboring N 
sites. A-mer correlations have been studied in various 
systems such as photons in dielectric waveguide [19, 20], 
acoustic waves [21] or electrons in crystalline lattice [22] . 
The off-diagonal disorder arises quite naturally, as the 
differences in the on-site energies change the tunneling 
amplitude values. As such, the off-diagonal disorder is 
correlated with the diagonal disorder and its relative am¬ 
plitude is small. On the contrary, in this paper we present 
a theoretical scheme allowing to control the off-diagonal 


disorder strength in an experiment (for the DRDM case 
see [23]). 

In Section II we describe a tight binding model with 
the binary disorder in both on-site energies and tunnel¬ 
ing strengths. Also, the methods of creating such a sys¬ 
tem using ultracold atoms in optical lattices is sketched. 
In Section III we describe a method employing transfer 
matrices to find resonant energies in the systems with 
binary correlated disorder, and use it to calculate posi¬ 
tions of all transparent states in the case of A-mers with 
no off-diagonal disorder. In Section IV we extend our 
analysis to include the diagonal - off-diagonal correlated 
disorder. We derive an equation for the resonant energies, 
discuss the asymptotic behavior and compare the results 
with numerical calculations of the localization length. 
Furthermore, in Section V we describe a dual random 
dimer model (being a 1-mer model with off-diagonal cor¬ 
related disorder). Again, we find [23] transparent modes, 
present an approximate expression for the localization 
length and compare the results with the numerical cal¬ 
culations. In the Appendix, a detailed description of the 
Floquet formalism used to derive effective Hamiltonian 
with off-diagonal disorder is presented. 

II. MODEL 

Consider a non-interacting tight-binding one dimen¬ 
sional Hamiltonian 

-ffo = X + h-c.)) , (1) 

i 

where ai(a|) is an operator of particle annihilation (cre¬ 
ation) at site i, rii = aj^a^ a particle number operator and 
Ci are on-site energies (we set h = 1 and use tunneling 
amplitude as an energy scale). The binary disorder oc¬ 
curs when {ci} are random numbers and take only two 
values: Eq or e^. Because of its simplicity, the binary dis¬ 
ordered models serve as a playground for studying the 
localization and transport properties. 

To create the binary disorder in an optical lattice one 
can use two species of atoms repulsively interacting with 
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each other. We allow only one fraction of atoms to move 
freely on the lattice, while the second one is trapped in 
the lattice sites (we call them frozen and denote with a 
/ superscript). The dynamics of the system is described 
by a Hamiltonian 

i? = ^ (Vbnfrii - +h.c.)) , (2) 


random dimer model (DRDM) where no two frozen par¬ 
ticles can appear on the adjacent sites. Such correlations 
can be created by using several lattices with different lat¬ 
tice constants (as described for the case of the DRDM in 
[24]). 

III. DELOCALIZED MODES FOR iV-MERS 


where we eliminated interactions between the mobile par¬ 
ticles (experimentally, one can switch them off by the 
optical or microwave Feshbach resonance [9] or use spin 
polarized fermions). As n-^ is fixed, we can treat Vo^f as 
an effective on-site energy e^. Furthermore, if the frozen 
particles are fermions or strongly repelling (hard-core) 

f 

bosons, then their occupation per a lattice site nl is ei¬ 
ther zero or one and hence the on-site energies take only 
two values Ei S {0, Vq}- 

The strength of the interspecies interaction Vb can be 
changed in an experiment by applying a proper magnetic 
field B, which is close to the value of the Feshbach res¬ 
onance. In our scheme, we choose the magnetic field to 
be time periodic B{t) = B{t + T), in such a way that the 
interactions are changing harmonically: 


A transfer matrix method is a very effective tool for 
the numerical computation of the localization length [25] . 
The formalism can be also used to find analytically the 
delocalized modes that are appearing in systems with 
the correlated binary disorder. A transfer matrix for a 
system described by (5) has a form: 



( 7 ) 


where E is an energy of a state. Iterating the procedure: 


( V’jAi 

V 


= T 


T 

1 . 


-ipi 

V'i-i 



( 8 ) 


Ci -t ei{t) = n{{Vo + Vi sm(ujt)). (3) 

As changing the magnetic field in a vicinity of the Fesh¬ 
bach resonance leads to rapid changes of the interaction 
strength, the relative values of Vb and Vl can be chosen 
at will in a broad range of values. In this way we obtain 
a time-dependent Hamiltonian: 

H{t) (ri{{Vo + Vl sm{ujt))ni - (ajoi+i -f h.c.)) . 

i 

( 4 ) 

If the modulation frequency w is bigger than other en¬ 
ergy scales in system {uj ^ 1) then, by the means of the 
Floquet theory, one can average out fast oscillating terms 
and find an effective Hamiltonian governing the long term 
dynamics of the system (see Appendix): 

■^e// = X! “ ti{a\ai+i E h.c.)^ , (5) 


where = n/Vb and ti = Jo — n{)j is a renor¬ 

malized hopping. For a binary case nf G {0,1}: 


f 0, if n/ = 0 
IVo, ifnf = l 



if n{ = n{+i 
if n{ ^ nf+i 


( 6 ) 


In our analysis, we consider systems with binary dis¬ 
ordered correlations in a form of Wmers where frozen 
particles always come in series of length N. All such 
blocks have the same internal structure - it is the dis¬ 
tribution of them which is random. Therefore, one can 
study transfer matrix of single block ranging from i + 1 
to i + N (i.e. T'l+^ = at •... • Tj+i) to get some insight 
in the transport properties of the whole system. For the 
case of an A^-mer with uniform tunnelings {t' = 1), such 
a matrix can be written explicitly: 


rpi-\-N 

-^2+1 


= T^{e) 


WN{e) -WN-i{e) \ .g, 

WN-l{e) -wn-2{£) ) ’ ^ 


where e = E — Vq and Wn{£) is a polynomial of the n-th 
order given by a recursive formula: 


Wnie) =-£Wn-i{e) - Wn- 2 {£), ( 10 ) 

where W-i{£) = 0 and w- 2 {£) = — 1. Solving the equa¬ 
tion (10) we can get an explicit expression: 



In particular, a transfer matrix for a system without 
disorder (for an arbitrary number of sites) reads 


where t' = Jo {Vi/oj) varies from 1 to about -0.4 (to the 
minimum of the Bessel function, min^, Jo(x)). 

Due to Anderson theory of localization in one dimen¬ 
sional system all states should be localized, unless dis¬ 
order have some correlations. In the binary disordered 
systems correlations can take a form of iV-mers where 
frozen particles always come in series of length N, or dual 


wm{E) -wm-i{E)\ 

^ \wm-i{E) -wm- 2 {E) j ' 


( 12 ) 


We can compare a free transfer matrix (12) with the N- 
mer transfer matrix (9) 

T^{E-Vo) =T^{E). (13) 
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If the condition (13) is met for some energy Er, then we 
can say that for a state with such an energy the A^-mer 
disordered system has the same transport properties as 
a free system. Hence, this state is delocalized. Such a 
procedure can be easily extended to any kind of correla¬ 
tions in a form of randomly distributed blocks with fixed 
internal structure, not necessarily fV-mers. 

A problem of finding transparent modes for binary 
disordered system with 7V-mer correlations can be fully 
solved by writing explicitly equations (13), using (10) and 
the fact that determinant of T^{e) is one (det(T'^(e)) = 
1). Then, it can be shown that the only nontrivial solu¬ 
tions (i.e. Vo 7 ^ 0) of the equation (13) exist if 


wn-i{E-Vo)=0, (14) 


resulting in (where we denote the matrix elements with 
vm, vm- 1 , vm -2 for the later convenience) 


r^e) = 


VM —VM-1 
VM-1 —Vm-2 


(18) 


where 


VM-2 = -^WNie), 

E 

VN -1 = -^WN{e) - Wjv-i(e), 

E^ 

VN = -^WN{e) + 2EwN-i{e) + t''^WN- 2 {e). 


(19) 


It is straightforward to check that a condition 


therefore a problem of finding resonant energies is re¬ 
duced to finding zeros of the polynomial wat-i- It imme¬ 
diately follows that: 

r^(A.-Ho)=±(j J), (15) 


f ^(e) = T^{E) (20) 

is not satisfied by energies given by (16). To fully solve 
the equation (20), we use the same procedure as in Sec¬ 
tion III. We assume that (19) fulfill a recursive formula: 


which means that (13) is trivially satisfied for M = 0. 
It turns our that a 7V-mer problem has — 1 resonant 
energies (this result, from a different approach, has been 
first presented in [26]): 

Er = Vo + 2cos{^i) for i e {1,..., A^ — 1}. (16) 

The additional condition for the appearance of extended 
states is that its energy has to lie in the energy band for 
a system without disorder {\Er\ < 2). 


IV. OFF-DIAGONAL CORRELATED 
DISORDER 

The properties of the system change dramatically if 
we add the off-diagonal disorder, as in (5). The off- 
diagonal disorder arises naturally in the model, as the 
differences in the on-site energies alter the hooping terms, 
but usually their values are small comparing to the di¬ 
agonal disorder. Nevertheless, in our approach we can 
independently manipulate the values of diagonal and off- 
diagonal disorder. For the V-mer, due to (6), only the 
edge tunneling amplitudes are changed: U = = t' 

(we assume here that V-mers are separated by at least 
one empty site). To obtain a transfer matrix of the N- 
mer, we need to write a transfer matrix for V -|- 2 sites 
(since the altered tunneling amplitudes t' appear in Ti, 
Ti_|_Ar and Ti_|_ 7 v+i we have to include two extra 
matrices): 

T N f — rpi-\-N-\-l _ rp T' T' T' T' T' 

^ i — i-\-N ■ • ■ -^i-\-2 

T^-2(£) 

(17) 


Vn = -EVn-1 - Vn-2- ( 21 ) 

It is a necessary but not sufficient condition for the reso¬ 
nance, as we do not know if there exist M G Z such that 
(20) is satisfied (in other words, we do not know if initial 
conditions of (21) meet that of (10)). 

From (21) we can get an equation for zeros of polyno¬ 
mial of N-th order in s: 

((1 - t'^)e - t'^Vo)wN-i{e) + (1 - t'^)wN- 2 {e) = 0. (22) 

First, we can analyze limiting cases: for F = 1 we restore 
the equation (14), so the approach is consistent with the 
one presented in Section III. In the case of t' —)■ ±oo 
we get wjv- 2 (e) = 0, so the resonances can appear for 
the same energies as for the {N — l)-mer with uniform 
tunnelings. The case of t' = 0 is trivial, as the chain 
is broken into pieces, and hence no transport can occur. 
Still, we can find that (22) comes down to wjv(e) = 0, 
so for t' —>■ 0 the resonant energies approach those of the 
{N + l)-mer. 

Solving (22) numerically gives always N solutions, for 
F —)■ 1 one of them diverges and the rest converges to 
N — 1 solutions for the V-mer with uniform tunnelings. 
If t' —)■ c» two edge solutions diverge and the rest of them 
go asymptotically to the solutions for the {N — l)-mer. 
An example for the 4-mer with Vq = 0.3 is presented 
in Fig. 1. Nevertheless, only for F = 1 we can say that 
states with Er move through the system in a way as there 
were no disorder. For t' ^ 1 situation is ambiguous, the 
equation (21) is satisfied, but we are unable to find M 
for which T^{e) — T^{E). However, for K —>■ oo, the 
difference asymptotically diminishes: 

min ||r^(£) - T^(A )||2 0. (23) 
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FIG. 1. (Color online) Blue solid lines: resonant energies for 
4-nier for Vo = 0.3 (solutions of (22)) in function of t'. Also 
band edges {E = ±2) and t' — 1 where solutions are given by 
(16) are indicated. 



FIG. 3. (Color online) Anderson localization length obtained 
from transfer matrix calculations for 4-mer model with Vo = 
0.4 and mean occupation of frozen particles 0.55. Red dashed 
line t' = 1, black solid t' = 0.9. 
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FIG. 2. (Color online) Asymptotic fall of matrix difference 
{rmuMKK |jT^(£) — r^^(i?)|| 2 ) to zero in function of K. Re¬ 
sults for 4-mer with Vb=0.5 and t' = 0.5 for lowest resonant 
energy Er = —1.1485. 


where ||.||2 is Frobenius matrix norm. Numerical results 
for the 4-mer with Vq = 0.5, t' = 0.5 (for Er = —1.1485) 
are shown in Fig. 2. 

In Fig. 3 we present numerically calculated Anderson 
localization lengths (using the standard transfer matrix 
method, see e.g. [25]) for the 4-mer and the mean occupa¬ 
tion 0.55, the strength of disorder is set to Vq = 0.4. The 
red dashed line corresponds to the system with uniform 
tunnelings f = 1 while the black solid line is obtained 
for F = 0.9. Resonances in case of = 1 are present 
at positions given by (16) while positions of resonances 
for t' — 0.9 agree with those obtained using (22). Values 
of localization length in the vicinity of Er for the case of 
t' — 0.9 shows that the resonance indeed exists. Localiza¬ 
tion length seems to diverge as increasing the resolution 
makes it bigger. 


V. DUAL RANDOM DIMER MODEL 


All the above analysis was valid for A-mers with N > 
2, the case of 1-mer is special. For the case of uniform 
tunnelings t' = 1 it does not have any delocalized modes 
- which is obvious as it is just the uncorrelated disorder. 
Even if we add an auxiliary condition that no two frozen 
particles can occupy adjacent sites (which is equivalent to 
a condition of separation of A-mers), situation does not 
change - it follows trivially from (13). However if t' 7 ^ 1, 
we get the so called dual random dimer model (DRDM). 
It has been extensively studied[23, 24, 27]. Such a model 
has the resonant energy 


Er 


Vq 

1 - t'2 ’ 


(24) 


for which: 


f\Er-VQ)=T^{Er) (25) 


is satisfied. As M is finite, this resonance is similar to 
resonances for A-mers with uniform tunnelings rather 
than for case with F yf 1. The resonance (24) is present 
as long as the resonant energy lies withing the band 
width < 2 which is equivalent to the condition: 

Vb < 2 1-F2|. 

The presence of the resonant mode results in the diver¬ 
gence of the AL length. In the limit of a small disorder 
(both diagonal Vq « 0 and off-diagonal t' « 1 ) it is pos¬ 
sible to obtain an analytical expression for the inverse 
localization length ]23] 


_ P (^0 + 2(1 - t'"^) cos{k)Y 
{1 + pY 8sin^(fcj 

X (1 _ ) , 

\ 1-b -b 2pcos(2fc) / 


( 26 ) 
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FIG. 4. (Color online) Anderson localization length obtained 
from transfer matrix calculations for 1-mer model with Vo = 
0.4 and mean occupation of frozen particles 0.55. Red dashed 
line t' = 1, black solid t' = 0.5 (DRDM). 

where we substituted E = 2cos(fc) for an eigenenergy 
E, and p = p/{l — p) where p is the mean occnpation 
nnmber of frozen particles. Indeed, the expression (26) 
vanishes for the resonant energy Er- 

VI. CONCLUSIONS 

In this work we have described a method of preparing 
the binary disorder in optical lattices nsing two atomic 
species. In snch systems, the off-diagonal disorder can 


be created using the fast periodic modulation of the in¬ 
terspecies interactions. Conseqnently, effective values of 
diagonal and off-diagonal disorder can be changed inde¬ 
pendently in a broad range. In order to observe single 
particle localization, we assnme that either onr mobile 
particles are spin polarized fermions or that we switch 
off interactions nsing the optical Feshbach resonance. 

We made a detailed analysis of the localization prop¬ 
erties in systems with the disorder with correlations in 
a form of 7V-mers, especially determining the resonant 
energies for which the nnobstructed transport through 
the system occurs. For that purpose, we employed the 
transfer matrix formalism. In this way, we restored a 
known resnlt in the case of no off-diagonal disorder [26] . 
Fnrthermore, we extended the analysis to the case with 
disorder in tnnnelings and derived an eqnation for the 
valnes of the resonant energies and analyze their asymp¬ 
totic values. 

We showed that, except the case of t' = 1, initial condi¬ 
tions of recurrence eqnation (10) are met only asymptot¬ 
ically. It is not clear whether these states are extended, 
but it seems that particles with the nearly resonant en¬ 
ergies always escape finite systems. 
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known as the Floquet states. Functions \un{t)) are 
T—periodic and e„ is called the quasienergy. Although 
are not eigenstates of H{t), it appears that \u{t)) 
are eigenstates of the Floquet Hamiltonian 

H{t) = H{t) — idt (A4) 


existing in the extended space of T-periodic functions. 
States in the extended space can be numbered using a 
new quantum number to S Z: \u^) = where 

|u™) is the eigenstate to the eigenenergy e™ = £„ J- com. 
It is straightforward to check that 

‘lO = (A5) 

so does not depend on the choice of to - adding 

w to £„ in physical space is equivalent to going into the 
next "Brillouin zone" for the quasienergies. Therefore, an 
eigenvalue problem for Hamiltonian (A4) can be reduced 
to a single "Brillouin zone", as long as we find its block- 
diagonal form in m-ordered basis. 

Let us perform a unitary transformation U: 


W = UHU\ 


U = exp 


^cos(ujt) 
—id - 

UJ 



(A6) 


Appendix A: Fast modulation of on-site energies 

In this section we derive an effective Hamiltonian for a 
cold atom realization of binary disorder with off-diagonal 
terms. The scheme we present allows to manipulate ex¬ 
perimentally relative values of diagonal and off-diagonal 
disorder. 

In Section H we consider the on-site energy variations 
due to disorder caused by a frozen fraction of particles. 
We apply periodic modulations of the on-site energies 
and obtain time dependent tight-binding Hamiltonian 

H{t) (n{{Vo Vi sm{ujt))ni - (al'oi+i J- h.c.)^ . 

(Al) 

As long as the driving frequency co is large (w 1) 
then we can expect that there exists an effective time- 
independent Hamiltonian iLe// which governs the long 
term dynamics. The explicit form of iLe// can be found 
very rigorously using the well established formalism of 
Floquet theory [6-8, 28-32]. 

The Floquet theorem, being a time analogue of the 
Bloch theorem, works for time periodic Hamiltonians 
H{t) = H{t -\- T). In this case, the time dependent 
Schrodinger equation 


idt\tpn{t)) = H {t)\'lpn{t)) 

(A2) 

has solutions of a form 


|^/’n(t)) = e“'®"‘|lt„(t)) 

(A3) 


U'{t) = ^ - (6**^"’+^ ’ alai+i J- h.c.)^ . 

(A7) 

Although after the transformation U a Hamiltonian ma¬ 
trix 

= (A8) 

is not block diagonal, it turns out that coupling between 
different Floquet blocks (to to') can be neglected in 
the high frequency limit (w 3> 1). Consequently, we can 
consider only one diagonal block which is responsible for 
the long term {t ^ 1/w) dynamics: 

1 

He« = f dtn'{t) = ^ - (tialm+i -f h.c.)^ , 

(A9) 

where 

ti = J7o ^ —(ci+i — (AlO) 

is the effective position dependent hopping and JTq is the 
zero-th order Bessel function. 

It is worth noting that U adds only local phases to sin¬ 
gle particles states and does not alter the density distri¬ 
bution. Therefore, localization properties of Hamiltoni¬ 
ans (Al) and (A9) are the same. Moreover, if the change 
of parameters is sufficiently slow, then it is usually possi¬ 
ble to adiabatically pass from eigenstates of one iLeff to 
another [33]. In particular, an eigenstate will follow the 
change of the disorder parameters during the experiment. 






